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Abstract 

Using a generalized cut vertex expansion we introduce the concept of an extended 
fracture function for the description of semi-inclusive deep inelastic processes in the 
target fragmentation region. Extended fracture functions are shown to obey a stan- 
dard homogeneous DGLAP-type equation which, upon integration over t, becomes 
the usual inhomogeneous evolution equation for ordinary fracture functions. 
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1 Introduction 



Within QCD, Wilson's Operator Product Expansion (OPE) [|TJ finds its most straightfor- 
ward application in lepton hadron Deep Inelastic Scattering (DIS). The cross section for 
this process is related to the time-ordered product of two currents, T[J(x)J(0)], which, 
at short distances, is expanded as 



The coefficient functions Ck{x) embody the short-distance structure of the amplitude, 
while Ok are local operators whose matrix elements describe long-distance effects. For 
DIS the OPE leads to an expansion of the forward Compton amplitude which, by use 
of dispersion relations ||, gives predictions for the moments of the DIS inclusive cross 
section. 

Unfortunately, for more general hard processes the straightforward use of the OPE 
technique fails. Two ways out of this impasse have been proposed. The first one || 
relies on a careful study of collinear and soft singularities order by order in perturbation 
theory and aims at proving that such singularities can all be lumped into some universal 
(i.e. hard-process independent) functions. This method has the advantage of being close 
to the physical parton picture for the hard process. As an alternative, an expansion in 
terms of cut vertices was proposed by Mueller 0] as a more direct extension of the OPE. 
Although more formal, the method of cut vertices has the advantage of giving directly 
an expansion for the absorptive part of the amplitude, i.e. for the cross section, in terms 
of the convolution of a cut vertex with a coefficient function. The cut vertex contains all 
the long-distance effects, and, as such, generalizes the matrix element of a local operator, 
while the coefficient function, as usual, is calculable in perturbation theory. The cut 
vertex method has been used in Ref. to deal with a variety of hard processes in QCD. 

Let us consider a deep inelastic semi-inclusive reaction in which a hadron in the final 
state is detected. When the transverse momentum of the hadron is of order Q the usual 
current fragmentation mechanism holds 0. In the region in which the transverse momen- 
tum of the produced hadron is much smaller than Q, the so-called target fragmentation 
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region, the simple description based on parton densities and fragmentation functions fails. 
For this reason, in order to describe the target fragmentation region as well, a new for- 
mulation of semi-inclusive hard processes in terms of fracture functions was proposed 
and later developed || [J. A fracture function M AA ,(x, z, Q 2 ) gives the probability of 
extracting a parton % with momentum fraction x from a hadron A while observing a final 
hadron A' in the target fragmentation region with longitudinal momentum fraction z. 
The possible relevance of such an object was already mentioned by Feynman |TU| before 
the advent of QCD. 

In this paper we will consider the case in which also the transverse momentum of the 
hadron A' (or equivalently the momentum transfer t = — {pa — Pa') 2 ) is observed. We 
will argue that, in the region i C Q 2 , an expansion in terms of cut vertices can be given 
and that a t dependent (or extended) fracture function^ can be defined in terms of a new 
kind of cut vertex. We will show, within a (0 3 )6 toy-model field theory, that extended 
factorization is just a consequence of collinear power counting [HI . 



The paper is organized as follows. In Sect. 2 we recall the definition of cut vertices 
and we obtain the expansion of the moments of the structure function in terms of matrix 
elements of local operators times coefficient functions by using cut vertices. In Sect. 3 the 
cut vertex approach is extended to deal with semi-inclusive processes and a generalized 
cut vertex expansion is given. In Sect. 4 we define extended fracture functions and their 
evolution equation is given and discussed. In Sect. 5 we discuss the results and draw our 
conclusions. 



2 Cut vertices 

In this section we will recall the definition of cut vertices in the simpler case of (0 3 )6 

theory. This toy model, despite its simpler structure, shares several important properties 

with QCD. It has a dimensionless, asymptotically free coupling constant and the diagrams 

with leading mass singularities have the same topology as in (light cone gauge) QCD. For 
1 Similar objects have been defined in the context of diffraction in Ref. |p"lf . 
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these reasons 



is an excellent theoretical framework for the study of factorization 



properties 



Consider the inclusive deep inelastic process 



p + J{q)^X 



off the current J = \$ 2 . We define as usual Q 2 and x as 



r\1 2 Q 2 

Q = -q x = — 

2pq 



F(j>,q) = ( f I fyr'"" < ,>\Jij))JiSy)\,>:>. 
In 



(2.1) 



Let us choose a frame in which p = (p + ,p_, 0) with p + ^> p_ and pq ~ p+q~- Given a 
vector k = (k + , k_, k) define k = (k+, 0, 0). The structure function, defined as 



(2.2) 



describes the interaction of the far off-shell current J(q) with an elementary quantum of 
momentum p through the discontinuity of the forward scattering amplitude (see Fig. 1). 

P q 




Figure 1: Deep inelastic structure function in 



The leading contribution to the structure function comes from the decomposition 
shown in Fig.2. Here, with the notations of Ref. 0, r is the hard part of the diagram, 
i.e. the one in which the large momentum flows, while A is the soft part. Decompositions 
with more than two legs connecting the hard to the soft part are suppressed by powers of 
1/Q 2 . 
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Figure 2: Relevant decomposition for the deep inelastic structure function in (</> 3 )6 



Such decomposition can be written in formulae as 



F(p,q) = J2 / V x (p,k)H T (k,q) 



afk 



(2.3) 



(27T) 



where V\(p,k) and H T (k,q) are the discontinuities of the long and short distance parts, 
respectively. 

Moreover, in order to pick up the leading contribution in eq. (|2.3| ), the momentum k 
which enters r is taken to be collinear to the external momentum p. Neglecting renor- 
malization, let us define for a given decomposition into a A and a r subdiagram 



Here v\(p 2 ,x) represents the contribution of A when the hard part is contracted to a 
point, while C T (x, Q 2 ) is the hard part in which one neglects the virtuality of the incoming 
momentum with respect to Q 2 . Since 




(2.4) 



and 



C T (x,Q 2 ) =H T (k 2 = 0,x,q 2 ). 



(2.5) 



Q 2 



(2.6) 



2p + q. 



and 



H T (k,q)=H T (0, 



Q 2 




(2.7) 



2k+q- 
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using the definition of k and eqs. (2.3)-(2.7) we can write 
r - d & k 

= E / ^(P 2 ' n ) ^" = / <V\ u) C(x/u, Q 2 ) ^. (2.8) 

The last integral defines the spacelike cut vertex v(p 2 , x) and the corresponding coefficient 
function C(x, Q 2 ). As usual a simpler factorized expression for the structure function is 
obtained by taking moments with respect to x. Defining the Mellin transform as 

U= I dxx^fix), (2.9) 
Jo 

we find immediately 

F a (p 2 ,Q 2 ) = v (T (p 2 )C (T (Q 2 ). (2.10) 
It was shown in Refs. [|15|, |16|] that the cut vertex represents the analytic continuation in 



the spin variable of a matrix element of operators of minimal twist. This correspondence 
has been confirmed up to two loops by direct calculation of the anomalous dimensions of 
cut vertices and leading twist operators [0. Hence, in the case of DIS, the factorized 



expression (p.10 ) can be identified with the one given by OPE 



F n (p 2 ,Q 2 ) = A n (p 2 )C n (Q 2 ) } (2.11) 

where A n (p 2 ) are now matrix elements of local operators. Thus, for integer values of a, 
the coefficient function which appears in (|2.10|) is the same as in Q2.ll ). This fact will 



be used in the next section where the evolution of the extended fracture function will be 
shown to be driven by the anomalous dimension of the same set of local operators. 

3 A cut vertex approach to semi-inclusive processes 

Let us consider now, still within (0 3 )6, a deep inelastic reaction in which a particle with 
momentum p 1 is inclusively observed in the final state, i.e. the process 

V + J{q) -> p' + X. 
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By using the same line of reasoning as for the inclusive case we may define a semi-inclusive 
structure function as (see Fig. 3) 




Figure 3: Deep inelastic semi-inclusive structure function in {4> z )q 



W(p,p',q) 



2vr 



[ d 6 x e iqx < p\J{x)\p'X >< X p'\J(0)\p > 

x 



(3.1) 



in terms of matrix elements of the current operator between the incoming hadron with 
momentum p and the outgoing hadron with momentum p' plus anything. 

When the observed particle has transverse momentum p'J_ of order Q 2 the cross section 
is dominated by the current fragmentation mechanism and can be written in the usual 
factorized way 



W(p,p',q)= / f A (x',Q 2 )a(x/x r ,z',Q 2 ) D AI (z/z',Q 2 ), 

J x' z 



(3.2) 



where 



pp' P- 
z = — ~ — 



(3.3) 



pq q_ 

In the language of cut vertices, eq. (|3.2j) is a convolution of a spacelike and a timelike cut 
vertex through a coefficient function fl 



W(p,p',q)= [ ^^v(p 2 ,x')C(x/x',z',Q 2 )v'(p>\z/z'). 
J X z 



(3.4) 
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By contrast, the limit t = — (p — p') 2 <C Q 2 is dominated by the target fragmentation 
mechanism and has not been considered in either approach^. Following the same steps as 
in the previous section, one can argue that, in the region t <C Q 2 , the leading contribution 
to the semi-inclusive cross section is given by the decomposition shown in Fig.4. 



Figure 4: Relevant decomposition for the semi-inclusive structure function in (4> 3 )e in the 
limit t <C Q 2 

Such a decomposition implies that an expansion similar to ( |2.8|) holds, in terms of a 
new function v(p,p',x) and a coefficient function C(x, Q 2 ) 

/A 
v(p,p',u)C(x/u,Q 2 )— (3.5) 
u 

where we have defined a new variable z as 

z = — ~ — (3.6) 
pq p + 

and a rescaled variable x = x/(l — z). The new function v(p,p', x) is given by 

v(p,p',x) = [ T{p,p',k)x8 fx ^—r) 7^-T6' ( 3 - 7 ) 

where T(p,p', k) is the discontinuity of a six-point amplitude in the channel (p — p' — k) 2 . 

The function v(p,p',x) is a new object that we will call a generalized cut vertex, which 

depends both on p and p' and embodies all the leading mass singularities of the cross 
2 In Ref. Q it has been shown at one loop that in the limit t — > a new collinear singularity appears 
in the semi-inclusive cross section which cannot be absorbed into parton densities and fragmentation 
functions, and so must be lumped into a new phenomcnological distribution, i.e. the fracture function. 
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section. By taking moments with respect to x as in eq. (|2.9| ), eq. Q3.5|) becomes 

W CT (p,p',q)=v a (p,p')C (T (Q 2 ) (3.8) 

that is a completely factorized expression analogous to fl2.10| ). 

We are now going to show that this expansion holds up to corrections suppressed by 
powers of 1/Q 2 - In order to do so we will use the method of infra-red power counting 



12, [L3J applied to our process. 




Figure 5: General reduced graphs which contribute to the semi-inclusive structure function 
in (0 3 ) 6 

In order to get insight into the large Q 2 limit of the semi-inclusive cross section, let 

us look at the singularities in the limit p 2 , p' 2 , t — > 0. The infra-red power counting 

technique can predict the strength of such singularities. Starting from a given diagram, 

its reduced form in the large Q 2 limit is constructed by simply contracting to a point all 

the lines whose momenta are not on shell. The general reduced diagrams in the large 

Q 2 limit for the process under study involve a jet subdiagram J, composed by on-shell 

lines collinear to the incoming particle, from which the detected particle emerges in the 

forward direction^, a hard subgraph H in which momenta of order Q circulate, which is 
3 In the large Q limit p and p 1 can be taken as parallel. 
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connected to the jet by an arbitrary number of collinear lines. Soft connections between 
J and H can be possibly collected into a soft blob S which is connected to the rest of the 
diagram by an arbitrary number of lines (see Fig.5). In (0 3 )6, by using power counting 
[13 1, we fmdQthat the leading contributions come from graphs with no soft lines and the 
minimum number of collinear lines connecting the hard to the jet subdiagram, as in Fig. 6. 




Figure 6: Leading contributions to the semi-inclusive structure function in (0 3 )6 

Any other diagram containing additional collinear lines between J and H is suppressed 
by powers of 1/Q 2 - It follows that W(p,p', q) is of the following form 

/d 6 k 
— T(p,p',k) H(k,q) + 0(1/Q 2 ). (3.9) 

It is now straightforward to show that eq. ( |3.9|) is equivalent to ( |3.5|) with the substi- 
tution H(0, x, q 2 ) = C{x, Q 2 ). Thus the expansion ( |3.5| ) corresponds to taking the leading 
part of the semi-inclusive cross section. 

4 This fact has been verified by an explicit one- loop calculation in Ref. ]l7|. 
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4 Extended fracture functions 



In the previous section we have given arguments for the validity of a generalized cut 
vertex expansion for the process p + J(q) — > p' + X in the region t <C Q 2 . We now want 
to investigate the consequences of such a result. 



The coefficient function which appears in ( |3.5[ ) is the same as that of ( |2.8| ) since it 
comes from the hard part of the graphs which is exactly the same as in DIS. So we can 
draw the important conclusion that the evolution of the coefficient function appearing in 
( |3.5| ) is directly related to the anomalous dimension of the leading twist local operator 
which drives the evolution of the DIS coefficient function. 

Despite the fact that the theoretical framework in which we have been working is the 



model field theory 



we expect that the main consequences expressed in (|3.5|) will 



remain valid also in a gauge theory such as QCD. The only further complications which are 
expected to arise are due to soft gluon lines connecting the hard to the jet subdiagrams. 
Unlike in (0 3 )6, in QCD these diagrams are not suppressed by power counting; hence the 
only way to get rid of such contributions is to show that they cancel out. This is the 
issue of a complete factorization proof, which is beyond the aim of this paper. As already 
argued in Ref. flTTl we do not expect that this complication will destroy factorization . In 
QCD, by using renormalization group, we have 



CUQ 2 ) = C l n (Q 2 /Qla s ) 



f 



Cl(l,a s (Q 2 )), 



(4.1) 



where Qq is the renormalization scale, ats = a s(Qo)i 7^ is the anomalous dimension 
matrix of the relevant operators, and an ordered exponential is to be understood. Thus 
we can write the analogue of eq.( |3.8|) in QCD as 



W n (z, t, Q 2 ) = M* n (z, t, Q 2 ) C;(l, a s (Q 2 )) 



(4.2) 



where, by following Ref. we have defined a t-dependent fracture function 

r a s(Q 2 ) , ffM' 



Mi(z,t,Q 2 ) = V:(z,t,Q 2 ) 



(4.3) 
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i/Qo 




Figure 7: Extended fracture function 

just in terms of a cut vertex V£(z,t,Ql) (see Fig. 7). 

Inverting the moments and expressing the extended fracture function in terms of the 
usual Bjorken variable x, one finds that M-\a>{ x -> z i Q 2 ) obeys the simple homogeneous 
evolution equation 

Q 2 ^M\ AI {^ z, t, Q 2 ) =Y,J\ ^ K v( u > a s(Q 2 )) M\ A ,(x/u, z, t, Q 2 ) (4.4) 

j 1 - z 

where Kij(u,a), defined as 

Kijiu, a) = — dn ^f(a) u~ n , (4.5) 

2« y i_ioo 

is the same DGLAP kernel which controls the evolution of the ordinary parton distribution 
functions. This result looks particularly appealing since it means that the evolution of 
the extended fracture function follows the usual perturbative behaviour. One may ask at 
this point how this result matches with the peculiar equation which drives the evolution 
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of ordinary fracture functions @: 



1 — z 



w^) n(xlu > Q)p;(u)D >' A \w^Y Q )- (46) 

The evolution equation for M A A ,(x, z, Q 2 ) contains two terms: a homogeneous term 
describing the non-perturbative production of a hadron coming from target fragmentation 
and an inhomogeneous term whose origin is the perturbative fragmentation due to initial 
state bremsstrahlung. As discussed in Ref. J7J, the separation between perturbative and 
non-perturbative fragmentation introduces an arbitrary scale but the fracture function 
itself M A A ,(x, z, Q 2 ) does not depend on it. 

It can be shown Jl8j that the evolution equation ( |4.6|) can be derived by an explicit 
calculation using eq. ( |4.4| ) together with Jet Calculus rules []l9fl . By defining in fact the 
ordinary fracture function as an integral over t up to a cut-off of order Q 2 , e.g. eQ 2 with 
e < 1: 

M AA ,(x,z,Q 2 ) = J dtM 3 AA ,(x,z,t,Q 2 ), (4.7) 

the inhomogeneous term in the evolution equation is obtained by taking into account the 
Q 2 dependence of the integration cut-off. 

Moreover, the interplay between the scales Q 2 and t has a sizeable effect in terms of a 
new class of perturbative corrections of the form logQ 2 /t. Such corrections are large and 
potentially dangerous in the region t -C Q 2 since they can ruin a reliable perturbative 
expansion. Those terms are naturally resummed into eq. (|4.3| ). For the extended fracture 
function these corrections do play an important role for understanding the dynamics of 
semi-inclusive processes in the kinematic region we have been considering here ITS] . 



5 Conclusions 



In this paper we have presented an extension of the cut vertex formalism to deal with semi- 
inclusive processes. We have shown that an extended fracture function can be defined in 
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terms of a new kind of cut vertex. Such an unintegrated, t-dependent fracture function 
obeys a DGLAP evolution equation with the same kernel as the one appearing in ordinary 
parton densities. This result is particularly evident in the cut vertex approach and agrees 



with the analysis made in a different context in Ref.[TT|. As a non-trivial check, the more 
complicated evolution equation for ordinary fracture functions can be obtained from 
the simpler one obeyed by extended fracture functions. 

Although our results have been obtained within the {4> 3 )e model, where a treatment 
based on mass singularities and infra-red power counting takes a much simpler form, 
we believe the extension to QCD to be harmless. Indeed, the picture that emerges for 
semi-inclusive hadronic processes in the target fragmentation region appears to be fully 
consistent with the intuitive QCD-corrected parton picture. 

In addition, a new class of logarithmic contributions appears at A 2 <C t <C Q 2 . Ex- 
tended fracture functions could give an additional tool to handle this interesting regime 
which certainly deserves further study |18| . 



We finally note that the relation between ordinary and extended fracture functions is 
quite analogous to the one between the experimentally defined/measured diffractive struc- 
ture functions F^ 3 \xjp, /3, Q 2 ) and F% (xjp, /3, Q 2 , t) [f20] which depend upon observing 
large rapidity gaps. Relating such quantities to fracture functions is another interesting 
issue for future investigation. 
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